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Abstract. A classical Calogero model in an external harmonic potential is known 
to be integrable for any number of particles. We consider here reductions which play 
a role of "soliton" solutions of the model. Wc obtain these solutions both for the 
model with finite number of particles and in a hydrodynamic limit. In the latter 
limit the model is described by hydrodynamic equations on continuous density and 
velocity fields. Soliton solutions in this case are finite dimensional reductions of the 
hydrodynamic model and describe the propagation of lumps of density and velocity in 
the nontrivial background. 
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1. Introduction 

The harmonic Calogero model (hCM) PQ [2] describes one-dimensional particles moving 
in the presence of an external harmonic potential and interacting through an inverse- 
square potential. The Hamiltonian of the model reads 

1 ^ 1 ^ .2 
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2 ^ 



JV 



— icjXj + ig 



(2) 

where Xj are coordinates of particles, pj are their canonic momenta, and g is the 
couphng constant. We took the mass of the particles to be unity. 

The model (classical and quantum) occupies an exceptional place in physics and 
mathematics and has been studied extensively [31 HI [3]. hCM similarly to other Calogero- 
Moser systems can be obtained by the Hamiltonian reduction of the system of non- 
interacting Hermitian matrices moving in external harmonic potential |3]. In this 
reduction the N coordinates of particles Xj appear as eigenvalues of simply evolving 
N X N matrix. The model is completely integrable and its solutions can be presented in 
terms of the eigenvalue problem for a finite dimensional matrix (see Sec. |2]for details). 

A remarkable fact is that the hydrodynamic limit — )■ oo of system ([T]) can be 
found exactly using the methods of collective field theory |6l 13 H] or using the methods of 
[HI Hn]. The hydrodynamic Hamiltonian can be written in terms of density and velocity 
fields as 



H 



dx p 



V 

~2 



91 
2 



(vrp^ - (9^ In y/pY 



iJ^X^ 



dx p - V — iux + ig{7Tp^ — dx In y/p) 



+ const, 



where p^ is Hilbert transform of p defined as a principal value integral 



P 



H 



dy 



piy) 



X 



(3) 
(4) 

(5) 



The density and velocity fields have a Poisson's bracket 

{p{x),v{y)} = 6'{x -y). (6) 

In this work we stress that the hydrodynamic form (3][6) can be used even for the finite 
number of particles / p{x) dx = N (see Sec. [6|. 

A goal of this paper is to find "soliton solutions" of the system (3][6). Corresponding 
soliton solutions of the Calogero model without an external potential {u = 0) are well 
known. A single solitons solution was found in [HI [12], and generalized to multi-soliton 
solutions in [lOj. 

Let us first explain what we mean by a soliton solution. Soliton is usually defined 
as "a pulse that maintains its shape while it travels at constant speed" . Obviously this 
definition does not make any sense in the presence of an external harmonic potential. 
Instead, we should talk about finite-dimensional reductions of an infinitely dimensional 
system (3]|6 ). Namely, if there is a solution of that system of the form p{x, t) = p{x; {zj}) 
(and v{x,t) = v{x;{zj})) so that the time dependence of p and v is reduced to M 
complex parameters Zj{t) {j = 1,2,...,M) with known time dependence, we call it 
an M-soliton solution. For example, in translationally invariant systems a one-soliton 
solution has a form p{x,t) = p{x — z{t)) with z{t) = vt which is consistent with the 
standard soliton definition. 
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The main result of the paper is the M-sohton solutions of (3]|6). It is presented 
in Sec. 7.3 The complex parameters Zj{t) of this multi-soliton solution in turn satisfy 
a "dual" Calogero model (25). Therefore, the complicated dynamics of an infinite- 



dimensional hCM (|3j) is reduced to an M-dimensional dynamics of complex Calogero 
system. We have to stress here that finding an explicit solution is still a non-trivial 
problem as one also has to relate initial conditions Zj{t = 0) of a dual Calogero system 
(25) with initial density and velocity profiles of ([3]). This is done implicitly in ( 72|[73 ). 
The derivations used in obtaining ( 72|[73 ) are very close to the ones used in [TO] . 

Remarkably, the finite dimensional reduction can also be performed in the finite- 
dimensional hCM ([!]) with particles. The evolution of ([T]) with finely tuned 
initial conditions can be described as a motion of few complex parameters Zj{t) 
{j = 1,2, ...,M) with M < N. This result is not published anywhere to the best 
of our knowledge and is another important result of this paper. 

The organization of the paper is the following. To introduce some notations and 
for the reader's convenience we start with a brief review of a solution of hCM ([T| in 
Sec. |2j We formulate a self-dual dynamical system which can be reduced to hCM in 
Sec. [3] A similar self-dual system has appeared before in Ref.[Tn] for a trigonometric 
Calogero model. We extend it to hCM. We show that this self-dual system allows for 
the reductions which correspond to soliton solutions of hCM. Several examples of such 
reductions are given in Sec. |4| In Sec. [5] we encode positions of hCM particles and their 
momenta by poles of meromorphic functions and derive equations for those functions 
using the approach of [HI [10]. We use these equations to rewrite the dynamics of hCM 
in hydrodynamic form in Sec. [6] and present soliton solutions in the hydrodynamic limit 
in Sec. [7| In concluding Sec. |8]we discuss possible generalizations of this work and some 
open questions. Some details of calculations are delegated to appendices. 



2. Solution of hCM with particles 



Here we briefly review the explicit solution of hCM (see Ref.[3] for review). In this 
solution the coordinates of Calogero particles Xj(t) can be found at any time as 
eigenvalues of a simple matrix Q{t). For the sake of brevity we do not discuss here 
neither a geometric meaning of the solution nor how this solution could be obtained 
(see |3]). Instead, we just introduce notations and give explicit formulas that we use in 
this work. 

Let us introduce the following N x N matrices: 



L = L ± iujX, where Lij = Pi5ij + (1 — 



Mi, = g 



(7) 
(8) 

(9) 
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These matrices depend on time through Xj{t) and Pjif) and satisfy important identities: 



[X,L] =ig{e®e^ - 1), 

Me = and e^M = 0. 

Here = (1,1,..., 1). 

It is straightforward to show that the equations of motion of hCM ([T]) 

9 O \ - 1 



Pi 



9x 



J fc=l (fc^j) 



are equivalent to the following matrix equations 

X + i[M,X] = L, 
L + i[M,L] = -u'^X 

or equivalently 

= -i [M, L^] ± iuL^ 

written in terms of L and M matrices usually referred to as a Lax pair. It immediately 

(17) 



(10) 
(11) 



(12) 
(13) 

(14) 
(15) 

(16) 



follows from (16) (see also Eq. B.l) that the following quantities 

h = Tr (L-L+)'= = Tr (L+L^)'^ 
are integrals of motion of hCM. Jo = Tr 1 = is the number of particles while 

h = Tr (L-L+) = 2nhCM (18) 

is the Hamiltonian ([T]) itself. The higher integrals of motion Ik, k = 2, 3, . . . are in 
involution, i.e. have vanishing Poisson's bracket with each other. The existence of a 
high number of conserved quantities is the result of integrability of hCM. 

One can also write the solution of hCM as an eigenvalue problem of a matrix 
which can be explicitly constructed from the initial positions and velocities of Calogero 
particles. Namely, the trajectories of particles are given by eigenvalues of the following 
matrix 



g(t) = X(0) cos(u;t) + w"^L(0) sin(wt). 



(19) 



Here the matrices X{0) and L[0) are constructed from initial conditions Xj{0), Pj{0) 
using definitions ( 7p ). 



3. Dual Calogero system and finite-dimensional reductions 



In this section we consider a complexified version of hCM ( p|l3| . We parametrize the 
complex momenta Pj by complex numbers zj so that the system ( 12p3 ) is rewritten 
as equations symmetric in Xj and Zj (see ( 20pl ) below). We refer to the obtained 
symmetric system as to a self-dual form of hCM. The self-dual form of hCM ( 20pl ) 
makes explicit the duality between particles Xj and excitations (parametrized by Zj) of 
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Calogero system. It is different from the action-coordinate duality explored previously 
in classical Calogero systems |13] . The self-dual system for the trigonometric Calogero- 
Sutherland model appeared previously in [in] (see Appendix C). It is transparent in the 
Hirota form (42) as a symmetry between tau-functions r~ and (see [9]). 

After introducing the self-dual form of hCM we consider different reductions of this 
system: reductions of the number of points Zj in a dual model and a real reduction {xj - 
real). Both of these reductions combined produce soliton solutions for an original hCM. 



3.1. Self-dual Calogero system 



Here we consider well and Zj cLS complex numbers. We introduce the 

following dynamic system: 

i . _ . = J2 + ^9Y1 (20) 



fc=i(fc^i) 

M 



n=l ^3 
N 



E 



Z'n 



--WY. ' (21) 

1 / / \ - ™ ■ 1 ^1 7 

m=l(m=f=n) J = i 

for Xj(t) with j = 1,2, ... ,N and Zn{t) with n = 1,2, . . . , M. We start with the case 
M = N. Let us note for future use that there is a connection between the motion of 
center of masses of points Xj and Zn obvious from ( 20|[21 ) 

N M 

y^jij - i^Xj) = ^(in - iujZn). (22) 

j=l n=l 

The system ( 20pT ) is Hamiltonian. It can be defined by its Hamiltonian given up to an 
additive constant —ugN{N + l)/4 by 



n 



hCM 



(23) 



+ 



9' 



N 



Xj Z]^ 



ug \- Xj + Zk 

Xj Z]^ 



E 

j,fc=l 



Y!ik=i Sjkdzk A dxj, where Sjk = ig{xj - Zk) 



j,k=i 

and by a symplectic form Q 
and corresponding Poisson's bracket {zk,Xj} = {S~^)kj. We notice that the system 
( [20|21|23[ ) is symmetric under simultaneous exchange Xj -H- zj and g — )■ —g. 

Equations ( 20|[2T ) are first order differential equations. The dynamics is fully 
defined by initial values of complex xj, Zj, i.e., by 2A^ complex numbers. 

Taking a time derivative of (20) (and similarly of (21 )) and using ( [20pl| we exclude 
first time derivatives. [|] As a result we obtain the decoupled systems of second order 

I After excluding first derivatives one has to reorganize products of fractions to exclude Zj from the 
first equation. For this purpose the following identity comes in handy ^^^^ + ^^^^ ^^^^ + = 0. 
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differential equations 



d 



2 dx 

2 dz.. 



N 
M 



2 UJ Xj, 



bj Z. 



31 



3 



1, 



AT 



M. 



(24) 



(25) 



The system (|24j) is a complex version of the system of equations of motion obtained 
from hCM ([T]), i.e., equivalent to ( 12p3 ). We refer to a system (25) as to the Calogero 
system dual to ( [24| ) or simply: the dual Calogero system. We outline the Lax formalism 
for this dual system and its correspondence to the one for the original system of Sec. |2] 
in 



Appendix A and Appendix B 



As soon as initial values of Xj and Xj are chosen, their evolution is defined by (24). 
Then the motion of complex points Zj is, on one hand, defined by the motion of Xj 
through ( 20||21 ) while on the other hand they evolve as Calogero system (25). This 
shows that one can think of the transformation Xj{t) — )■ Zj{t) given by (20) as of the 
Backlund transformation from one solution of (24) to the other. We do not explore the 
connection of our results with Backlund transformations further in this work. 



3.2. Reduction of number of particles in a dual system 

A remarkable fact is that the derivation of ( 24p5 ) from ( 20pT ) also holds if M 7^ 
(we are interested here in M < A^) and one can still think of (25) as of a dual system 
for (24) consisting of smaller number M < A^ of particles. The difference with M = N 
is that in the case M < N one can not generically solve (20) to find Zn for an arbitrary 
choice of Xj, Xj. Some fine tuning of initial values of Xj, Xj is necessary. Instead, one can 
specify A^ complex points xj and M complex points Zn and then find xj from (20). Then 
by (21) the motion of A^ points Xj is reduced to a motion of M < A^ complex points 
Zn governed by a dual Calogero system (25) having fewer degrees of freedom than the 
original system (24). We refer to this reduction as to a dimensional reduction or to an 
M-soliton reduction of (24). 

The soliton reduction can also be understood as a limit in which some coordinates 
of dual particles go to infinity. Indeed, let us consider the self-dual system ( 20||2l ) with 
M = N. We choose initial positions xj arbitrarily and initial positions Zn so that the 
latter are divided into two groups. The coordinates z„ {n = 1,2,..., M) are arbitrary 
while the coordinates Zk {k = M + 1, . . . , N) are very far away from the origin so that 
for k > M: \zk\^ \xj \ for any j and \zk\ ^ \zn\ for any n < M. We are interested in the 
limit Zfc — )■ 00 for k > M. One can see that in this limit only M coordinates Zn, n < M 
enter the equations for xj as it is written in (20) with M < N. The equations for Zj are 
divided in this limit into M equations for i„ with n < M (see (21)) and to completely 
decoupled system of A^ — M points Zk [k = M + 1, . . . N). The latter system is not 
important for us while the system ( 20||21 ) with M < N gives an M-soliton reduction 
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as the dynamics of Zn with n < M is given by (25) having less degrees of freedom than 



(24). 



3.3. Real reduction 



So far we considered Xj as complex numbers. It is clear, however, from (24) that once 



initial values of Xj and Xj are chosen to be real they stay real at later times, even though 
Zj are moving in a complex plane. Let us specify some arbitrary real values of Xj{t = 0) 



and Xj{t = 0). For M = N one can generically solve an algebraic system (20) (A^ 
algebraic equations with M = N unknowns zj) and find corresponding initial complex 
Zj(t = 0) and then using (21) initial Zj(t = 0). This procedure defines a real reduction 



of the complex system ( 20pT ). We can think of ( [20pl| ) as an alternative way to write 
the system ( 12p3 ) or equivalently (24) understanding that initial complex values of Zj 
and not arbitrary but constrained by reality of Xj and Xj. 



We notice here that a true symmetry between (24) and (25) exists only for complex 



variables Xj and for M 



N. Imposing reality conditions on x 
the symmetry between (real) Xj and (complex) z„. 



explicitly breaks 



4. Soliton solutions of hCM with N particles 

Now we consider the case when both real and soliton reductions are applied 
simultaneously. In this case one can take real and imaginary parts of complex equations 



(20) and write the following real equations 

1 



UJX-i 



9 



N 
M 



M 

E 

n=l 



+ 



■ 9 



71 = 1 



Xj Zfi 



(26) 



(27) 



If M complex positions Zn are given at any time one can find both real positions Xj and 
corresponding real momenta pj. The data Xj,pj are not independent but "tuned", i.e., 
related by ( 26p7 ) through the values of M complex parameters z„ (2M real parameters). 



Equations (26) have an electrostatic interpretation. Indeed, (26) can be obtained 



as extrema conditions for the following function 

N 

E 



N 2 
UXj 



^ N M 

- ^ In - Xfcl + - ^ ^ In \xj - Zn\^. 

7=1 n=l 



E = - l^\Ti\x, - xu\ + -l^l^ln \x, - z^\\ (28) 

j=i ^ j<k 

This function coincides with an "electrostatic energy" of particles with unit charges 
interacting through a logarithmic potential (2d Coulomb potential). The particles are 
restricted to move along a straight line (a real axis) and are in the presence of 2M 
external charges —1/2 placed at Zn.,Zn and an external harmonic potential. We notice 



here that the solution of (26) is not necessarily a minimum of (|28j). Soliton solutions 
correspond to any extremum (maximum, minimum or saddle point) of (28). It is 
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important to stress that here and in the following we choose the signs > and 



(7 > which guarantees that the harmonic potential in ( 28 ) is confining. 



4-1. Background 

As an ultimate case of M-soliton reduction we consider M = which gives a static 
solution. Indeed, ( 26p7 ) in the limit — > 00 for all n becomes pj = for all j and 



coordinate of particles in equilibrium are defined by (26) as: 

^ 1 
= 9 J] • 



(29) 

fc=l(fc^j) ■ 

It is well known that a solution of this system of algebraic equations is given by the 
roots of A^-th Hermite polynomial (Stiltjes formula [IHIIS])- Namely, 



h 



11 



HN{h,)=0, J = 1,2, 



N. 



(30) 



4-2. One soliton solution 

Consider the case M = 1. Equations ( [26p7l ) give 



UJX-i 



Pj 



9 



N 

E 

.9 I 1 



+ 



X 



1 



(31) 



(32) 



The equations (31) can be viewed as a new generalization to the Stieltjes problem 



(29) (see Refs. [m [161 Ej). To the best of our knowledge this generalization to the 



Stieltjes problem has not been studied and exact solutions of (31) are not known. One 
can think of (31) as of definition of some polynomials Pn{x,zi) = Yiji^ 



such 



that PN{xj,zi) = for j = 1,2,...N. In the limit — )■ 00 we have Pn{x,zi) — >■ 
H]\j{x\Ju] / g). We make some progress in describing these solutions in the limit ^ 1 
in Seclll 

1 takes an especially simple form 

(33) 



The equation (25) in the case M 



and can be easily solved 

zx{t) 



i.e., the trajectory of Zi is an ellipse in the complex plane. Using (22) for M 
obtain the parameters of this ellipse 

sin cut 



zrit) = z^{Q)e'^' + 



[P(0) - iujX{0)] 



(34) 
1 we 

(35) 



where -2i(0) is the initial position of zi in the complex plane and X = 
P = Ylf=iPj center of mass and the total momentum of the system at t = 0. 

Both A:(0) and P(0) are in turn defined by zi{0) through ( [3Tp2| . 
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Let us consider for simplicity a particular initial value ^i(O) = ih with 6 > 0. Then 



the solution of (31) gives X = 0|§J The equation of the ellipse in this case is 

z{t) = ibcos{ut) - (6 - P(0)/w) sin(wt), (36) 
where we find from ([32]) 

^(0) = -^E:;:2f^<0- (37) 



3 



The inequality means that a = b — P{0)/uj > b so that the major semiaxis a of the ellipse 
is always along the real axis. In the limit b oo, -P(O) ~ and major and minor 
semiaxes are ^ ~ ^ + ^ and b respectively. The eccentricity of the ellipse goes to zero 
(ellipse becomes a circle) as 6 — t- oo. In the opposite limit 6 — )■ we have -P(O) ~ — f 
giving ct ~ In this limit the ellipse has a large eccentricity with the major semiaxis 
a ~ — oo as the minor semiaxis 6 — )■ 0. 

Let us now fix some large value of the major semiaxis a by taking -2i(0) = a > 0. It 



is clear from this analysis that there are two different solutions of (31) corresponding to 



large and small values of minor semiaxis of the ellipse. These two solutions correspond 



to two different extrema of electrostatic energy (28). For one of them all N particles 
( "cloud" ) are located around the origin, far from the external negative "soliton" charge 
placed at a. For the other extremum the cloud around the origin consists of — 1 
particles. One more particle is far away from the cloud, close to the external charge. 
The former solution corresponds to the large minor semiaxis b ^ a, while the latter 
corresponds to b ^ g / a. If we decrease a two corresponding values of b approach 
each other and become equal to some "critical" value b = be- At this value the major 



semiaxis a has a minimum value Oc. For a < Oc there are no real solutions of (31 ). Later 



in Sec. [t] we will show that in the limit of large this minimum occurs at b^. ~ N^^^ 
and corresponds to a minimal major semiaxis Oc - ~ A^-l/^ where R = ^2gN/ U IS 
the radius of the "cloud" of particles. A world-line diagram of a typical single-soliton 
solution for b < b^ is shown in Figure [TJ In this regime the soliton solution looks like a 
Newton's cradle. The soliton is essentially a single particle when its position is outside 
of the "cloud". This particle transfers its momentum all the way through the system 
with the other particle being kicked out from the other side of the system. Due to 
the interactions (in contrast to the actual Newton's cradle) the particle is dressed by 
other particles when inside of the cloud. This picture was qualitatively described by 
Polychronakos jl8j . 



5. Particles as poles of meromorphic functions 

In this section we, following an approach of [U [10] consider particles of hCM as poles of 
meromorphic functions and derive dynamic equations satisfied by those functions. 



We do not know how to prove this statement. Equations (31 1 have a symmetry Xj — > —xj and 



numerical solutions suggest that this symmetry is unbroken resulting in X = 0. 
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Figure 1. World-line diagram for a single-soliton solution of Harmonic Calogero Model 
is shown for = 40 and the value 6 « 0.84 < 6^ ~ 1.67. Lines are the world lines 
of individual Calogero particles. There are no crossings of world lines. However, the 
lump of density corresponding to a soliton travels all the way through the system. It 
becomes an isolated particle outside of the "cloud" of other particles. 



We start by introducing two meromorphic functions u^{x) of a complex variable x 



N 



U [X] 



M 



+ iux, 



n=l 



X 



(38) 



(39) 



These functions are completely defined by their poles xj and z„ which move as Calogero 
particles ( 24p5 ). The function u~{x) is defined by its poles - the coordinates of hCM 
particles xj. The function u~^{x) is defined by the coordinates of the dual model Zn or 
alternatively by its values at xj given by 

^ 1 

u~^{xj) = pj — iuxj + ig —. (40) 



OC rj Ob If' 



Conditions (40) are equivalent to (20). Notice that the r.h.s. of (40) appears in the 
factorized form of the hCM Hamiltonian (|2|. 

Having defined u^{x) by ( [38p9| ) we can rewrite the system (|20[2T| as a single 
equation 

,2 



Ut 



U 

~2 



I- [u 



u 



oj'^x'^ 







(41) 



2 V '^2 

with u = ^- . Indeed, assuming the form ( 38p9 ) and taking the residues of (41) 



at points Xj, Zn we reproduce ( 20pl ) respectively. The equation (41) is a version of a 
bidirectional Benjamin-Ono equation |I0] modified for hCM. A key advantage of (41) 
is that the number of particles does not enter the equation explicitly and, therefore. 
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this form is well-suited for taking hydrodynamic limit N ^ oo. Before discussing this 

(42) 



limit in Sec. |6] we also give a bilinear Hirota form of (41 ) 

^2 



9 
2' 



OJxDr — 



T 



0, 



where are given by: 



u 
u 



—ig dx In r~ 



(43) 
(44) 



and, e.g., Dtf ■ g denotes Hirota derivative. We note here that up to trivial time- 
dependent factors tau-functions are given by 



N 



T (x) 

r+(x) 



~ ^j) — det(x — Q), 

M 

]^(x - Zn) = det(x - Q), 



(45) 



(46) 



n=l 



where the N x N matrix Q is given by (19) and Q is the corresponding dual M x M 
matrix (A.6). The self-duality of hCM is expressed then as an obvious symmetry of (42) 
under the exchange of tau-functions r~ r+, — )■ —g. 



6. Equations of motion in hydrodynamic form and hydrodynamic limit 



Here we rewrite equations of motion for hCM in a hydrodynamic form for finite and 
then consider the hydrodynamic limit of those equations, i.e., the limit of infinitely many 
particles N ^ oo. We start with equations for u^{x) and with corresponding analyticity 
and reality conditions and then present the equations of motion in hydrodynamic form, 
i.e., written for particle density and velocity fields. We again follow the approach of 

m- 

Let us start by rewriting hCM ([!]) in terms of fields 
identical to 

1 



u 



One can show that ([T| is 



n 



hCM 



/dz 

® dz 
4vr^ / 



u 



(47) 



u u 



4- -2 4.0 - 

+ U U H U OzU 



I 2 2- 
+ OJ Z U 



by using the definition (38) and the property (40). The contour of integration in (47) 
goes around the real poles Xj of u~{z) counter-clockwise and does not encircle any of 
complex poles Zn of u^{z). The equations of motion (41) is equivalent to ( T2p!3 ). 

The poles of u~{z) are real and one can parameterize the real analytic function 
iu^ {z) by a real function of a real variable p{x) - a particle density field. We introduce 

N 

p{x) = Y,6{x-x,) (48) 
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and rewrite ( 38 ) Cauchy transform of p{x): 

r> + 00 

u [z] = lujz — ig I dx 



p{x) 



Z — X 



(49) 



where z is a complex number not coinciding with any of poles of u~{z) (e.g., Im {z) ^ 0). 

The field W^iz) is discontinuous on a real axis with the discontinuity related to the 
density of particles. More precisely 

u~ {x ± iO) = ^'ngp + i{ngp^ + lux) (50) 

with the discontinuity 

u~ {x + iO) — u~ {x ~ iO) = —2'Kgp{x). (51) 

Using ( p8p9| as well as ( [20pl] ) and ( [50l|5l| after some calculations we obtain that on 
the real axis 

p{x)u'^{x) = —ig Ittp{x)p^ (x) — 2^xp] + iu!xp{x) + ''^^Xj6{x — Xj). (52) 
We identify the last term of the r.h.s as a momentum density of the system 

N 

p{x)v{x) = '^^Xj5{x — Xj), (53) 

where v{x) is the velocity field. We divide ( [52| by p{x) and obtain 

u~^{x) = V — ig [tcp^ — dx log ^/p) + icux. (54) 

Equations ( 50p4 ) give the relation between fields u^{x) and microscopic density and 
velocity fields. We notice here that these relations are exact even in the case of finite 
number of particles A^. The density and velocity fields have a conventional Poisson's 
bracket ([6]). Substituting ( 50p4 ) into the Hamiltonian (47) we arrive to the Hamiltonian 
of hCM in a hydrodynamic form ([3]). 

Hamilton equations following from (3]|6) are the Euler and the continuity equations 
for density and velocity fields: 

Pt + dxipv) = 0, (55) 



vt + dx{ — + w{p) + 



2 2 
where a chemical potential w{p) is given by: 

1 



(56) 



(57) 



We remark here that although the equations in this section are written in hydrodynamic 
form, they are still valid for a finite number of particles and are equivalent to the 
corresponding equations for hCM. In the case of finite N the density and velocity fields 
are singular functions given by their microscopic definitions ( |48p3| . All expressions 
involving these fields and their products should be properly regularized as it is explained 
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above. The key point of the regularization is to use the definitions ( 38p9 ) of as 
meromorphic functions. 

Let us now go to a hydrodynamic hmit. This simply means that from now on we 
treat p{x) and v{x) as continuous (even smooth) fields forgetting the discrete nature 
of hCM particles. Note, that the information about the total number of particles is 
still preserved in the relation J pdx = N and one should do some rescaling of fields 
when going to the large limit (see Sec. 7.2). Having specified an initial configuration 
p{x,t = 0), v{x,t = 0) one can in principle solve ( 55p6 ) and find density and velocity 
fields at all times. An interesting class of solutions (multi-soliton solutions) of ( 55|[56 ) 
is realized for a fine-tuned initial configurations of fields. As the number of particles 
A^ and the number of poles M of the field u~^{x) are independent parameters, one can 
take a hydrodynamic limit N ^ oo keeping M finite and fixed. As a result one obtains 
solutions in which the dynamics of continuous fields p{x,t) and v{x,t) is reduced to a 
motion of M points Zn{t) in a complex plane. This is a finite-dimensional reduction 
of an infinitely dimensional hydrodynamic system. We refer to this reduction as to an 
M-soliton solution. 

In the next section we consider several examples of soliton solutions in the large A^ 
limit. 



7. Soliton solutions of hCM in hydrodynamic limit 

7.1. Background 

Let us find the configuration p{x) and v{x) with given J p{x) dx = N that minimizes 
the energy We rewrite (|3]) in a manifestly positive form 



H 



dx p 



r^2 ^ 2 
y + 2 i^9P^ - gdx log ^/p - ^x) 



+ const 



(58) 



/ 



1 2 

dxp-\u^[x)\ + const. 



where we used (54) to obtain the last line. It is easy see that the minimal energy 
condition is 



M+(x) = 0, 



or writing it separately for real and imaginary parts and using (54): 

gij^p^ - log Tp) - = 0, 
v{x) = 0. 



(59) 

(60) 
(61) 



Eq. (60) is the hydrodynamic form of the equation (29). It describes the distribution 
of zeros of Hermite polynomials HN{xy^ujJg). In the limit A^ — )■ oo we think of p{x) 
and v{x) as of continuous fields. In this limit the solution of (60) is given by a Wigner's 
semi-circle law 



Tig 



x^ 



R 



2gN 



(62) 
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Eq. (60) also appears in the context of Random Matrix Theory (see, for example, Refs. 
[T9| 120]). We notice here that both the density at the origin po(0) = p 



1 



and 

9 

the radius of the cloud of particles R are proportional to N^^"^. The main correction to 



(62) in the next to leading order in comes from the fact that the largest zero of 
Hn{x\Ju} / g) is not R but is given asymptotically by Xmax ~ -R — where the 

constant 71 ~ 1.8557 ... is related to zeros of Airy functions. It is also notable that the 
distance between neighbor roots goes as Xn+i — Xn ^ ^ N'^/"^ close to the origin 
and xtv — a^Tv-i ~ R~^ near the boundary of the cloud. 



7.2. One-soliton solution 

The one-soliton solution is given by 



u'^ix) 



X 



with zi{t) satisfying (25) for M = 1 or (33). Using (54) we rewrite (63) as 



V — lUX 



(63) 



(64) 



X — Zi 

This relation allows one, in principle, to find density and velocity fields from the position 
zi at any moment of time. The soliton parameter zi{t) is moving in a complex plane 



along the ellipse (35). Therefore, (64) gives a 1-dimensional reduction of an infinite 



dimensional Calogero system in hydrodynamic limit defined by (3]|6). Eq. (64) is a 
hydrodynamic analogue of (20) with M = 1. 



Taking real and imaginary parts of (64) we obtain hydrodynamic counterparts of 



( |31p2| 



^(vrp^ - 9^ log ^) 



UJX 



X 



+ — 1 , (65) 

(66) 

X — Zl X — Z\ ) 

It is remarkable that the velocity field of a one-soliton solution is given explicitly by a 



2 



1 



simple expression (66). The equation (65) defines, albeit implicitly, the density field for 



a one-soliton solution. Comparing ( 65p6 ) with the corresponding background equations 
( |60p[| ) we see that the fields for a one-soliton solution are obtained by perturbing the 
background configurations by terms ~ 1/zi. In particular, in the limit 2:1 —t- 00 we go 
back to the equilibrium configuration ( 60|[6T ) . In the large N limit the term dx log ^Jp 
and the right hand side of (65) are both suppresed by 1/A^ with respect to other two 



terms. We also notice here that the solution of (33) for z\{€) is given by (35), where 



dx pf , 



X 



dx px. 



(67) 



are the total momentum and the center of mass of the system. Of course, finding P{0) 
and X{0) from zi{0) using ( [65p6| is still a non-trivial problem. 



CONTENTS 



16 



In the limit a; — )■ 0, — )■ oo and p = po(0) 



2uiN 

a 



const the equation 



(65) gives rise to Lorenzian shaped sohtons in agreement with solitons obtained by 
Polychronakos [H] and Andric et. al. [12] for a model without harmonic potential and 
with the background density p. 

As the exact solution of ( 3lp5 ) is not available we briefly discuss the solution in 
the limit of large A^ in the next to leading order in 1/A^. 

Rescaling variables (p, x, — )■ {p, x, Zi,v)y/N in one can easily see that 

the right hand sides of ( 65p6 ) are of the order of 1/A^. Therefore, in the leading order 
in A^ one has density and velocity given by ( [6lp2| ). 

The correction to (60) consists of two parts: the correction to the background 
solution without solition and to the correction caused by the presence of the soliton, i.e. 
by the right hand side of (65). Here we are interested only in the latter. 

First, let us assume that the solution of (65) is given by a smooth function p{x). 
Then we have: 

1 yi 



p{x,t) 
v{x, t) 



7T {x — Xi] 



[x - xiY + y{ 



2 ' 



(68) 
(69) 



where we denoted zi(t) = xi{t)+iyi(t). The solution ( 68|[69 ) describes a lump of density 
of the changing width ~ yi{t) located at the moving point Xi{t). The point Zi{t) moves 
according to (35). Let us start with zi{0) = ib. Using ( 68|[69 ) and (67) we find the 
parameters X(0) = and P(0) ^ -^+ujb- ujy/WT¥. The major semiaxis of the 



ellipse is given (see Sec-4.2) by 
a = b-uj'^P{0) 



9 
2ub 



(70) 



As a function of time zi{t) = ib cos ut — a sin cut. We can identify several interesting 
limits corresponding to different values of b. 



Large: b ^ R \/N. In this case we have a «i 6 + B?/{2b). The trajectory of Zi is 
close to a circle of a very large radius b. The width b of the soliton is bigger than the 
size of the cloud. In this case there is no pronounced lump of the density. Instead the 
whole cloud oscillates slightly around the origin. 

Intermediate I: N^^^ ^ 6 ^ \/N. For b R we obtain from 

The major semiaxis a has a minimum at be ^ {gR/2uy/^ ~ N'^/^^g/2uj. The width 
of the soliton is 6 at t = 0. It is much smaller than the size of the system and one can 
see a very well pronounced lump of density while soliton travels through the system. 
The width of the soliton somewhat changes in time but remains much larger than an 
interparticle distance inside the cloud. Therefore, the continuous approximation is still 
valid in this regime at all times. The soliton in this regime is a well-pronounced lump 
of density which oscillates inside the cloud of particles. 



70)a^i?+^ + fi. 
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Figure 2. Time evolution [ti < t2 < < of a one-soliton solution in hCM 
is shown schematically in the large N limit. The figure corresponds to the regime 
Intermediate II. As soliton moves to the left its width decreases and becomes of the 
order of interparticle distance at some x (shown by the dashed line). After this point the 
continuum approximation is not valid and soliton is represented by the delta-function 
(shown by an arrow) . This figure is a continuous analogue of Figure 



Intermediate II: N"^^"^ <^ b <^ N^/^ . This is, probably, the most interesting regime. 
For an initial configuration the width of the soliton h is still much bigger than an 
interparticle distance A^~^/^ in the middle of the harmonic trap. Therefore, we still can 
use a continuous approximation and the value ~ -R+ + 1^' However, as a function 
of time yi{t) decreases and at some point becomes of the order of an interparticle 
distance at the point Xi{t). [j]] Starting from this time we cannot use the continuous 
approximation. Instead, we assume that the density can be divided into a delta function 
corresponding to a single particle plus a continuous background with — 1 particles. 
In the limit when yi is much smaller than an interparticle distance we have simply 



in this regime a boundary particle is kicked out of the cloud and travels outside of the 
cloud for some fraction of the period of the motion. 

II A maximal interparticle distance is at the edge of the cloud and is of the order N^^/^ . [14] 
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Small: b <^ N ^^'^ . In this limit the continuous approximation is invahd already at 
t = and we consider an isolated particle at the origin with other A^ — 1 particles forming 



a continuous cloud (71 ) for all times. The value of -P(O) is given by a microscopic formula 
(37) which is dominated in this case by a particle at the origin P(0) ~ —g/b. It gives 
a ~ 6 + g/{(jjb). The density evolution is given by a particle moving in the semicircle 



background ( 7Tp2 ) 



7.3. Multi-soliton solution 

Here we briefly list equations describing M-soliton solutions of ( 55||56 ). The density and 



velocity fields are completely defined by complex coordinates Zn through (39) which can 



be rewritten using (54) as 
V — iujx = ig 



E 



X 



+ ig{7Tp -d^\ogy/p). 



(72) 



An initial configuration of M complex numbers Zn defines initial velocity and density 



fields through ( 72 ) . After density and velocity fields are found one can easily determine 



initial z- velocities using a hydrodynamic analogue of (21) which has a form 

M ^ 

Zn-iujZn = ig ^ — + 7ig{p + ip^). 



(73) 



After initial velocities i;„ are obtained, the dynamics of Zn is defined by (25) so that z„ 
can be found as eigenvalues of the matrix Q (A.6). 



The dynamical problem of finding a multi-soliton solution of ( 55p6 ) is reduced to 



finding the density p{x,{zn}) from (72). The latter is still complicated, albeit time- 
independent, integral equation. 



8. Conclusion 



In this work we used a self-dual formulation ( 20pT ) of a harmonic Calogero system 
(hCM) Q to find an M-soliton reduction of an hCM with A^ particles M < N. Soliton 
solutions can be obtained by fine tuning the initial conditions Xj, Xj for Calogero 
particles. We found a hydrodynamic formulation of this reduction and then took a 
hydrodynamic limit N ^ oo keeping M finite. As a result we obtained an M-soliton 
solution of an infinitely-dimensional hydrodynamic system. 

The derivations of this article are based on similar derivations of [10] made for 
Calogero- Sutherland model. We emphasized in this work that the soliton reductions 
are possible even for finite number of particles while Ref. [10] considers only finite 
soliton solutions of an infinite-dimensional translationally invariant model. We gave 
the generalizations of finite dimensional reduction to the cases of trigonometric 
and rational Calogero-Moser systems in Appendix C and Appendix D| respectively. 
The generalization of finite-dimensional reductions to the elliptic Calogero models is 
also rather straightforward and will be given elsewhere. We also expect that the 
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generalizations of our results to more general Calogero-Moser systems related to Lie 
algebras are possible. 

In this work we did not discuss the meaning of the presented soliton reductions 
of hCM neither within the projection method of solving hCM (see [Sj [23]) nor within 
inverse scattering formalism |2T] . It is interesting to find the corresponding descriptions 
of both self-dual formulation of hCM and of the soliton reductions. 

The self-duality of hCM given by ( 20pl ) and used in this work is different from the 
known dualities of Calogero models [22], [13] . It would be interesting to have a precise 
relation between those dualities as well as the connection with the known bispectral 
property of Calogero-Moser systems [211 [2S1 12S] • 

hCM and many other Calogero models remain integrable after quantization. 
Moreover, many results obtained for classical Calogero models have direct analogues 
for their quantum counterparts. In particular, the pole ansatz can be extended to the 
quantum case [9] . The classical soliton solutions of Calogero models correspond to quasi- 
particle excitations of the corresponding quantum models [5j. It would be interesting 
to give quantum analogues of the results presented here. 
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Appendix A. Lax formalism for a dual Calogero system 



In this Appendix we describe Lax matrices for a dual Calogero system. The formalism 
is almost identical to the one presented in Sec. |2] for an original hCM. In addition 
we introduce an intertwining matrix F relating corresponding matrices between dual 
systems. 

Let us define matrices dual to ( 7|8p ) as 



L ± iuZ, where L^n 
1 



9 



M 

E 

l=l(ljLm) 



ig 



(A.l) 
(A.2) 

(A.3) 



Similarly to (To|TT) we have: 

[Z, L] = ig{eM ® - 1a/), (A.4) 
Mcm = and e^M = 0. (A.5) 
Here e^^ = (1, 1, . . . , 1) is a row vector made out of M ones and 1m is a unit M x M 



matrix. To avoid confusion we will denote the vector e from (11) by here 
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It is obvious that other formulas of Sec. [2] can also be written in terms of dual 
variables and matrices. For example, similarly to Xj{t) the values of Zn{t) at any time 
can be found as eigenvalues of a Q matrix 

Q{t) = Z{0) cos{ut) + uj'^L{0) sm{ujt). (A.6) 

The dual variables Zj are related to original variables Xj through ( 20pT ) and the 
natural question is how the corresponding matrices and, in particular, integrals of motion 
for the dual system are related to the ones for an original system. Here we relate 
the matrices (7]|8) to ( A.1[|A.2 ). Then in Appendix B we find the relations between 
corresponding integrals of motion. 

In the following we consider matrices ( [7p|A.l|A.2[ ) as functions of parameters Xj , Zn 
only, with time derivatives expressed in terms of these parameters using ( 20pl ). We 
introduce one more "intertwining" rectangular matrix F of the size N x M 

Fin = (A.7) 

Xi Zn 

The matrix F depends on both direct and dual variables and provides a connection 
between dual systems as we will see below. It is straightforward to show that the 
following identity holds for both upper and lower signs 



L^F = FL^ - ujg{l ± l)eN ® eij, 
where multiplication is the matrix multiplication. We also find the identities 



L ejy = Fe 



F 



(A.8) 

(A.9) 
(A.IO) 

which are equivalent to equations of motion ( 20pT ). 

We conclude this Appendix by stating that the matrix F obeys a simple matrix 
evolution equation 

F = -icoF -iMF + iFM, (A.ll) 

while L"^ satisfy equations fully analogous to (16) 

± r . . n 

-i M,L^ ±iujL^. (A.12) 



A derivation of (A.ll ) is based on ( 20pT ) and is rather straightforward albeit somewhat 
cumbersome. 



Appendix B. Integrals of motion 



It follows from ( 16 ) that 



(B.l) 



and, therefore, (17) are integrals of motion. In fact, the time evolution (B.l) does not 
change eigenvalues of L~ and describes an isospectral deformation of this matrix. 
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Similar conclusion can be derived for matrices L and L^L using (A. 12). Here we 



relate the integrals of motion (17) to analogous expressions for the dual system. 



Let us start with an easily verifiable identity 

[L"^,L~] = 2iu[Z,L] = —2ug{eM 



-1m). 



We proceed as 



FL+L- 



- 2ujgF{eM 
L'(^L^F + 2ujgeN'S) 
L-L+F + 2ujgF, 



- Im) 



-A/ 



j - 2ujgF{eM ® - 1a/) 



(B.2) 



where we used (A. 8) and (A. 9). If / is an eigenvector of L^L^ and A is a corresponding 
eigenvalue i.e., L^L^ f = A/, it follows from (B.2) that L^L^ has an eigenvalue X + 2ujg 
with the corresponding eigenvector Ff. We conclude that M eigenvalues of L^L^ are 
identical (after the shift by 2ujg) to the M eigenvalues of L^L~ . We show below that 
the remaining N — M eigenvalues of L^L^ are constants given by (B.6). Therefore, 
integrals of motion of the original and dual hCM are simply related. 

We start with the relation between the integrals of motion of dual systems for the 
case M = N. In this case the matrix F is square and invertible (we assume that Zj ^ Xk 
for any j,k = 1, . . . A^). Then one can find the matrices for the dual system from 



(A. 8) etc. In particular, (B.2) can be written as 

= F(L+L- - 2uglN)F-\ (B.3) 

One immediately concludes that integrals of motion of dual systems are connected by a 
very simple relation 

4 = Tr{L-L+f = Tr{L^L- - 2w^l^)^ (B.4) 

To consider the case M < N we exploit the fact that the dimensional reduction to 
M-soliton solution can be obtained by taking some of Zj to infinity as it is described in 
Sec. 3.2 We divide zj into two groups. We keep zi, . . . ,zm finite and take ZM+k = Zk for 
A; = 1, . . . {N-M) to infinity. We take this limit for the matrix L+L- = {L-f + 2iuZL~ 
and leave only non-vanishing matrix elements. We use the fact that all Xj are chosen 
to be finite. The matrix obtained in the limit has a block-triangular form and its 
eigenvalues are given by the eigenvalues of {L^L~)m reduced to the size M x M and to 
the eigenvalues of the (A^ — M) x (A^ — M) matrix 2ugB defined as 



N-M 



N-M 



Zi 



It is easy to showj^ that eigenvalues of B are 1, 2, 3, 
N-M eigenvalues of L~ are trivial and given by 



5, 



i'3 J - 

Zi Zq 



(B.5) 



A^ — M. Therefore, the first 



\s = 2ugs, 



s = 0,l,2,...,{N-M-r 



(B.6) 



% The matrix B is triangular in the basis of f^''\ A: = 0, 1, . . . , (iV - M - 1) defined by (/('')); = (z,)''. 
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The remaining M eigenvalues are not trivial and coincide with those of the M x M 
matrix of the dual model L^L~ shifted by 2ujg. This fact illustrates the meaning of M- 
dimensional reduction for integrals of motion. In particular, for the background solution 
M = all iV eigenvalues of L~L^ are given by 2ugs, s = 0, 1, 2, . . . {N — 1). The latter 
result is known and can be obtained directly from the properties of Hermite polynomials 
(see eqs. 10a,b of Ref. [27]). 



Appendix C. Solitons as finite dimensional reductions of N-particle 
Sutherland Model 



Here, for the sake of completeness we give a self-dual form of the Calogero- Sutherland 
model (trigonometric case of Calogero model) as it appeared in [TOj. Then we give 
formulas for soliton reductions. 

Calogero-Sutherland Model describes particles on a circle interacting with inverse 
sine-squared (chord-distance) interactions 

N , ^ N n 



H 



E 



9 



sm 



(C.l) 



where L is the circumference of the circle. Positions and momenta of particles on a 

27TiXj 

circle can be characterized by Wj = and pj = —i{L/27r)wj/wj, where < Xj < L. 



The self-dual form of the Sutherland Model analogous to ( |20pl| is: 

N 



2tt 

T 



2tx 

T 




,k=l J 



k=l ^ 
M 

^ Uj 

k=l ^ 



Uk 



Uk 



(C.2) 



(C.3) 



for M = N. Here the "positions of solitons" are labeled by complex numbers uj with 
\uj\ 7^ 1. The finite dimensional reduction of the Sutherland model, i.e. M-soliton 
solutions are given by ( C.2[C.3 ) with M < N. 

Taking real and imaginary parts of (C.2) we obtain the following relations between 
soliton positions and positions and momenta of particles: 



N 



E Wj + Wk _^ i 
in ■ — 11), 9 ^-^ 



Wj 
= 1 ^ 



Pj 



■Kg_ 
'2L 



M 



M 



Wj + Uk 



+ 



Wj + 



k=l 



Wj - Uk Wj - 



E 

k=l 



Wj + Uk 
Wj - Uk 



Wj + ^ 
Wj — 



(C.4) 



(C.5) 



w. 



for particles on a circle. The static solution is obtained for 
M = 0. It is easy to check that up to translation it is given by pj = 0, Xj = jL/N (or 



Here we used that Wj 



Wj 



e N 



')• 
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Appendix D. Soliton reduction of Calogero model (rational case) 



Here we discuss how the soliton reduction can be implemented for the rational Calogero- 
Moser system or Calogero model (CM). This model is given by Hamiltonian ([T]) with 
a; = 0. It can be written in a self-dual form by taking w = in ( 20pT ). Then an M- 
soliton reduction can be obtained by taking M < in ( 20pT ). Although this reduction 
is well defined for a complexified system, applying it to the original real Calogero model 



we run into the following difficulty. The real equations (26) do not have solutions for 
M < A if = 0. It is easy to understand from the electrostatic interpretation. Indeed, 
it is not possible to keep A repelling charges within some finite interval on a line with the 
using the negative charge M < A in the absence of an additional harmonic potential 
(see (28) with w = 0). We show here how to overcome this difficulty and obtain an 



M-soliton reduction for CM. 

Let us consider the following change of variables 



(D.l) 
(D.2) 



cos (cut) 
ut' = tan(a;t). 

It is known that this transformation "removes harmonic potential" [3]. Namely, if Xj{t) is 
a solution of hCM, the transformed functions x'j{t') defined by ( D.1|D.2 ) give a solution 
of the Calogero model. 

It is clear that an M-soliton reduction of hCM gives through the change of variables 
( D.1[D.2 ) a corresponding reduction of CM. 

Let us apply the change of variables ( D.1[D.2 ) to the self-dual form of hCM ( 20pT ). 
We obtain 



lOOXj 



iujf 



tUJZ' 



1 + iut' 



M 



1 



n=l 



N 



where we also changed Zn — )■ z'^^ similarly to (D.l). We consider ( D.3|D.4 ) as a modified 
or deformed self-dual form of CM. uj is just a parameter of the deformation (there is 
no time scale u in CM). At the value a; = equations ( D.3[D.4 ) give an unmodified 
self-dual form of CM. At a; = there are no real solutions Xj{t) for M < A as it was 
explained above. However, for a; 7^ one obtains all soliton reductions corresponding to 
the ones for hCM. The obtained soliton solutions will have an explicit time dependence 
additional to the time-dependence of parameters Zn{t). 

Before giving an example of the reduction we stress that excluding z^s from 
( D.3[D.4 ) one arrives to the system of second order differential equations for CM. The 
parameter u does not enter these equations. Similarly, excluding Xj's one finds that the 
parameters Zn{t) form a dual CM, that is also move according to CM equations. 

Let us consider the simplest example of soliton solutions for CM. Namely, we 



CONTENTS 



24 



consider M = 0-soliton reduction of rational CM corresponding to a static (background) 



solution of hCM (30). This solution is mapped to 



xUt') 



(D.5) 



This equation gives 0-dimensional reduction of CM system. It is easy to check that, 
indeed, (D.5) solves (D.3) for M = 0. The parameter u enters the initial conditions 
[f = 0) of (D.5) and defines the time scale. The limit a; — )■ is singular and does 



not correspond to a physical solution of CM. In this Appendix we showed that soliton 
reduction of the the rational Calogero model can be implemented via mapping soliton 
solutions of hCM onto solutions of CM using ( D.1|D.2 ). The same kind of mapping can 
be done between two hCM with different frequencies which will result in new soliton 
solutions that will have additional explicit time dependence. 
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